Strong bounds on Onsager coefficients and efficiency 
for three terminal thermoelectric transport in a magnetic field 
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For thermoelectric transport in the presence of a magnetic field that breaks time-reversal 
symmetry, a strong bound on the Onsager coefficients is derived within a general set-up using three 
terminals. Asymmetric Onsager coefficients lead to a maximum efficiency substantially smaller 
than the Carnot efficiency reaching only in the limit of strong asymmetry. Related bounds are 

derived for efficiency at maximum power, which can become larger than the Curzon-Ahlborn value 
?7c/2, and for a cooling device. Our approach reveals that in the presence of reversible currents the 
standard analysis based on the positivity of entropy production is incomplete without considering 
the role of current conservation explicitly. 
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The Onsager reciprocity theorem derived more than 
eighty years ago is arguably the most important con- 
straint for the phenomenological framework of linear irre- 
versible thermodynamics [ll] . In its original form it states 
that the linear kinetic coefficient iife(B), which relates 
the affinity Tk to the flux Ji is equal to the reciprocal 
coefficient B) provided both i and k refer to quan- 

tities with even signature under time-reversal The 
constant magnetic field B breaks the time-reversal sym- 
metry. A generalization which includes variables with 
odd signature was given by Casimir in 1948 [sj. Ever 
since, the reciprocal relations have turned out to be ex- 
tremely useful for a plethora of applications. A promi- 
nent one is the unified theory of thermoelectricity which 
allows to treat the various thermoelectric effects such as 
the Peltier, the Seebeck or the Thomson effect on an 
equal footing thus revealing their interdependencies 0, Hi- 
lt is well known that properly designed energy filters 
allow a substantial enhancement of the efficiency of ther- 
moelectric devices fsl]. Linke and co-workers 0,0] showed 
that, in principle, even Carnot efficiency is attainable. 
Despite these promising theoretical results, the actual 
performance of present devices is notoriously much more 
modest, which keeps the search for better thermoelectric 
materials a very active and increasingly important re- 
search field as reviewed in [J H | . Recently, in an intrigu- 
ing paper, Bcncnti et al. pointed out that the presence 
of a magnetic field could in principle enhance the perfor- 
mance of a thermoelectric device [12!|. In fact, their quite 
general thermodynamic analysis invoking only the On- 
sager reciprocity relations and the positivity of entropy 
production suggested the possibility of a device deliv- 
ering finite power while still reaching Carnot efficiency. 
Such a spectacular option deserves both further scrutiny 
and a search for a microscopic realization. In this Let- 
ter, we show that unitarity of the scattering matrix as 
a general physical principle imposes a strong restriction 
on the Onsager coefficients that lead to a significant re- 
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duction of the attainable efficiency of any thermoelec- 
tric device within the broad and well-established class of 
ll-[l6|. 

transport can conveniently be dis- 
cussed within the model sketched in FigH] Two particle 
reservoirs of respective temperatures Tl > ?r and chem- 
ical potentials /il < are connected by a conductor C, 
which allows for the exchange of heat and particles. Con- 
sequently, as soon as the steady state is reached, constant 
heat and particle currents, Jq and Jp, flow between the 
left and the right reservoirs. Depending on the sign of 
these currents, the machine works either as a power gen- 
erator or a refrigerator. We now assume that both the 
temperature difference AT = Tl — Tr > and the chem- 
ical potential difference = /iL — /^r, < are small 
compared to the respective reference values, which we 
choose to be r = Tr and /i = /iR. In this linear response 




FIG. 1. Sketch of a thermoelectric device in the presence of 
a magnetic field B. The entire dashed box represents the 
conductor C, which is connected to two reservoirs. For the 
special case of the three terminal model, the conductor con- 
sists essentially of a scattering region and an additional probe 
terminal. 
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regime, the currents Jp and Jq are related to the affinities 
Fp = A^/T and = AT/T^ via the phenomenological 
equations 



Jq — LqpJ~ ~^ LqqJ' q. 



(1) 



The constant rate of entropy production accompanying 
this transport process reads 

S — J-pJ p-\-J-qJq — LppJ-p -\-LqqJ-q -\'{Lpq-\-Lqp ):FpFq. (2) 

Clearly, the features of this machine are completely deter- 
mined by the kinetic coefficients Lij — p,q), which 
are subject to two fundamental constraints. First, the 
second law requires 5 > 0, which is equivalent to the 
conditions 

Lpp, Lqq > and LppLqq - {Lpq + Lqpf/A > 0. (3) 

Second, Onsager's theorem implies the symmetry rela- 
tion 



Lpq{B) = Lqp{-B 



(4) 



Here, we have reintroduced the magnetic field B which 
was suppressed so far in order to keep the notation slim. 

In the light of the Onsagcr reciprocity relation (U), 
the expression ^ for the total entropy production rate 
suggests a natural splitting of each current J,; into a re- 
versible and an irreversible part defined by 



jicv _ ^ij 



L 



'3^ 



(5) 

Obviously, J^^^ vanishes for B = 0, i.e. for unbroken 
time-reversal symmetry. However, once the latter is bro- 
ken, although not contributing to 5', the reversible cur- 
rents can, in principle, become arbitrarily large, since be- 
sides (jSj no further general relations between Lij(B) and 
Lji(B) are known. Such unconstraint reversible currents 
could ultimately give rise to the possibility of dissipation- 
less transport noticed by Benenti et al. [12| as mentioned 
above. 

For progress investigating possible further constraints 
the phenomenological set-up discussed so far must be 
made more specific. A paradigmatic model for the con- 
ductor C is shown inside the dashed box in Fig. [TJ It 
consists of a central scattering region, with a constant 
magnetic field B, and a third electronic reservoir (ter- 
minal) P, whose temperature and chemical potential are 
chosen such that no net exchange of particles and heat 
with the rest of the system occurs. The scattering re- 
gion is connected to all of the reservoirs via perfect, 
one-dimensional, infinitely long leads. In order to keep 
the model as simple as possible, we assume moderately 
damped @, non-interacting electrons, which means that 
no inelastic scattering events take place outside the reser- 
voirs and the electrons are transferred coherently between 



them. This set-up constitutes the quite general class of 
three-terminal models on which we base our subsequent 
analysis. We emphasize that the additional terminal P 
plays a crucial role since it is well known that in a purely 
coherent two terminal set-up the off-diagonal Onsager 
coefficients must be even functions of the magnetic field 
17] and hence there are no reversible currents. Such 
probe terminals, whose temperature and chemical poten- 
tial are chosen self-consistently, were originally proposed 
by Blittiker Ts'l and have become a common tool to simu- 
late inelastic events in an otherwise conservative system. 

Since we now have to deal with three reservoirs, we 
need accordingly four affinities J-^ = {^a — fJ-)/T and 
J^^ = {Ta - T)/T^ {A = L,P), where we stiU use the 
temperature and chemical potential of the right reservoir 
as reference values /x and T. For convenience, we collect 
the affinities in two vectors = [F^^F^). Analo- 



gously we define the current vectors J = {eJ^,Jq), 



I 'Jq 

with J^ and J^ the particle and heat currents fiowing 
out of reservoir A and e the electronic unit charge. This 
vector notation allows us to write the phenomenological 
equations in the rather compact form 



with 



^PL '^PP ' 



Here, the 1^'ab (^j B — L, P) are 2 x 2-matrices of kinetic 
coefficients. By virtue of the additional constraint = 
0, we can moreover eliminate and thereby reduce ([5]) 
to a system of two equations 



with an effective matrix of kinetic coefficients 



^LP- 



'PL 



(7) 



(8) 



relevant for the net currents from L to R. 

Explicit expressions for the 2x2 block matrices 
of kinetic coefficients showing up in arise from the 
multi-terminal Landaucr formula [19, |20| 
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AB 



where 
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[SAB-\SABiE,B)\'), (9) 



F{E) 



E- n 
kBT 



(10) 



is the negative derivative of the Fermi function, ks de- 
notes Boltzmann's constant and h Planck's constant. 
Sab{E, B) 7^ Sba{E, B) are the matrix elements of the 
3x3 scattering matrix S(E, B) that describes the passage 
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of electrons with energy E through the scattering region. 
Current conservation requires that S(i?, B) is unitary and 
the time reversal invariance of unitary dynamics implies 
the symmetry S{E,B) = §(£;,-B)* We emphasize 
that in the presence of a magnetic field, this symmetry 
permits B) to be non symmetric. Consequently, de- 
spite the fact that the are still symmetric, the re- 
duced matrix ([8]) will generically acquire a non- vanishing 
asymmetric part as desired. 

We now derive a constraint on the asymmetry of the 
reduced matrix L. First, we define the Hermitian matri- 
ces 



(L'+L")-f-iV3(L'-L'*) 



(11) 



and 



IK EE (L + L*) + iV3{h - L*) = (^Pf ^"^^ , (12) 

\^pq ^qqj 

where L' is the full 4x4 Onsager-matrix introduced in 
while L is the reduced matrix Second, as a conse- 
quence of the unitarity of the scattering matrix B), 
K' can be shown to be positive semidefinite on |2l| . 
Third, since for any z e 
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it then follows that the reduced matrix K must be pos- 
itive semidefinite on C^. Therefore the matrix elements 
of K have to obey the inequalities 



Kpp.Kqq > 



and 



KppKqq 



KpqK;^ > 0. 



FIG. 2. Bounds on benchmarks in units of i]c as functions 

(13) of the asymmetry parameter x. The upper diagram shows 
^max, the lower one ?7(Pmax)*. The solid curves represent 
the bounds following from relation pSI) . In the limit x — 
±oo both functions asymptotically approach the value 1/4, 
which is shown by the dash-dotted lines. For comparison, the 
bounds obtained by Benenti et al. [12(1 solely from the second 
law Q have been included as dashed curves. The dotted line 

(14) in the lower panel indicates the Curzon-Ahlborn limit 1/2. 



We note that the Kij — p,q) are rather complicated 
functions of the matrix elements of the full Onsager- 
matrix L', which would make it a quite challenging task 
to obtain directly from ([5]) and If we insert the 
definition (1121) of the Ka in terms of the Li, , we imme- 



diately get the relation 



LppLqq 



LpqLqp 



r2 



Lip > 0, 



(15) 



which constitutes our first main result. We emphasize 
that ()15|) provides a much stronger constraint than (|3]) 
since the former can be rewritten as 



The new bound (|T5|) has profound consequences for 
the performance of the model as a thermoelectric heat 
engine. Following the lines of Benenti et al. 12|, we 

introduce the dimensionless parameters 



y = LpqLqp/Det L and 



Lpq/Lq 



h{x) < 4y < if a: < 0, 
< 4y < h{x) if x>0, 



(17) 



Expressed in terms of x and y, the inequality (1151) reads 

(18) 



LppLqq 



{Lpq+Lqp)y4>3{Lpq~Lqpf/4. (16) 



Indeed, we recover ^ only iiLpq is equal to Lqp. As soon 
as the kinetic coefficients contain finite asymmetric parts, 
the left hand side of (fTS)) must be strictly larger than zero. 
In other words, the reversible currents associated with 
the asymmetric part of L come at the price of a stronger 
lower bound on the Onsager coefficients and hence on the 
entropy production rate ^ than the bare second law (|3]) 
requires. 



where h{x) = 4,x/{x — 1)^ . The most important bench- 
marks for the performance of a heat engine - namely the 
maximum efficiency Tymax and the efficiency at maximum 
power 7?(Pmax) " admit the analytic expressions [l^ 



.W = .cx^^^and,(P_) 



Vc- 



xy 



2y 



(19) 



Here, we have denoted by 77c = 1 - Tr/Tl w TT^ the 
Carnot efficiency, which is the absolute upper bound for 
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FIG. 3. Bounds on the coefficient of performance of a re- 
frigerator e in units of rjQ as a function of the asymmetry 
parameter x. The dashed hne follows from bare second law, 
the solid line from the stronger relation (|15[l . 



the attainable efficiency following from the second law. 
Both benchmarks become maximal for 4y = h{x). The 
resulting bounds, r;*^^^ and ?7*(Pniax), which are our sec- 
ond main result, are plotted in Fig. [2l It shows how 
the maximum efficiency decays rapidly as the asymme- 
try parameter x deviates from its symmetric value I. In 
the limit x — > ±00, maximum efficiency and efficiency at 
maximum power both approach rjc/^- This result im- 
plies essentially that from the perspective of maximally 
attainable efficiency the thermodynamic cost of the re- 
versible currents is larger than the benefit they bring. 
For efficiency at maximum power, the situation is some- 
what different. If x is only slightly larger than I, the 
Curzon-Ahlborn-limit r]cA — ?7c/2 j22l426l |. reached for 
X = 1, can be overcome in a small range of x values with 
a maximum of 47^(7/7 at x ~ 4/3. This result shows that 
despite the strong bounds on ?7max it may in principle be 
possible to improve the performance of the machine by 
breaking the time-reversal symmetry slightly. 

Finally, we discuss the consequence of the bound (fTS]) 
for the model as a refrigerator. In this case the most im- 
portant benchmark is the coefficient of performance e 



-J; 



defined as the ratio of the heat current 
extracted from the cold reservoir and the absorbed power. 
In terms of the dimensionless parameters ([T7|) e reads 12 1 



X vV+T+l' 



(20) 



where 77^ = Tji/{Ti, - Tr) w l/TJ"^ is the efficiency 
of an ideal refrigerator. Like for r^max and ?7(Pinax), the 
maximum s* of e is attained for Ay = h{x). The resulting 
bound is plotted in Fig. [H Again we observe that the 
efficiency deteriorates rapidly as x deviates from 1. 

In conclusion, we have investigated the thermoelec- 
tric transport properties of the most general version of 
the paradigmatic three terminal model in the presence 
of broken time reversal symmetry. We have derived a 



strong constraint on the linear transport coefficients that 
can be obtained neither from Onsager's principle of mi- 
croreversibility nor from thermodynamic arguments in- 
voking only the second law. This constraint implies 
strong bounds on both maximum efficiency and efficiency 
at maximum power if this device is operated as a ther- 
moelectric heat engine or refrigerator. 

We emphasize that our results rely solely on the uni- 
tarity of the scattering matrix. Therefore our new re- 
lation (jlSp applies to any specific model which can be 
formulated within the three terminal set-up, including 
especially the ones based on quantum dots jla ] and bal- 
listic micro junctions [l5| . which have been proposed in 
this context recently. It might even apply to a model 
introduced by Entin-Wohlman and Aharony IJ], which 
includes a phonon bath locally interacting with electrons 
transferred between the reservoirs. Furthermore it is 
noteworthy that by following exactly the same lines the 
constraint (fT5|) can be shown to apply to the three ter- 
minal railway switch transport model [l3|. This model 
may be regarded as a classical analogue of the quan- 
tum three terminal model discussed here. Indeed, our 
considerations are not restricted to the quantum realm, 
since ultimately the unitarity of the scattering matrix is 
nothing but the manifestation of the law of current con- 
servation, which should be considered the fundamental 
principle underlying our bound. 

Despite this generality, the question whether a similar 
bound can be found for more complex model classes, 
including for example a larger number of probe terminals 
or models with genuinely interacting electrons, remains 
open and should constitute an important subject for 
further investigations. 
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Here, we prove that the matrix 

K' = (L' + L'*) + iV3{h' - L'*) 



(1) 



as introduced in Eq. (11) of the main text is positive 
semidefinite on C*. To this end, we first notice that by 
virtue of the multi-terminal Landauer formula [cf. Eq. 
(9) and Eq. (10) of the main text] 



(5AB-\SAB{E,m- 



(2) 



the matrix elements K'l^-^ = {K')im (Z, to = 1, . . . , 4) can 
be written in the form 



with 



)' 111' 



and the Hermitian matrix 



(4) 



2(1-|^LLP) 
SlP? + \SPL? - lVi{\SLP? -\SPL?] 



\SlP? + \SPL? + iVi{\SLP? - \SPL?) 

2(l~|5pp|2) 



(5) 
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m 



X 
J3 



From (m onwards we notationally suppress the explicit 
dependence of the v; , V and the elements of the scatter- 
ing matrix on energy E and magnetic field B. For any 
z = (zi, Z2, z-i: ZiY G we have 



z^K'z = zfzmK'i^ 

l,m=l 

= "^fjE F{E) zt.\\j V z,^.^| (6) 



of the determinant yields 



dE F{E)y^Yy 



with 



E 



(7) 



The last line of ([6]) reveals that if V is positive semidefi- 
nite, the same holds for K'. 

In order to show that V is indeed positive semidefinite, 
we first observe that both diagonal entries of V are obvi- 
ously non-negative. It therefore suffices to show that the 
determinant of V is non-negative. The direct calculation 



Det V 



il-\SLL\')il-\Spp\') 

-\Slp\^-\Spl\^ + \S, 



5LPI \Spl\^. (8) 
Now, we distinguish two cases. First, we assume 
I^LPp > \Spl\^ and rewrite dH]) as 

^ = {\SLP\' + \SLHmSLP\' + \Snp\') 



\Slp\^-\Spl\^ + \Slp\'\Spl\^ 

\SpLf{\SLpf -\SpLf)>0. 



\Slp\^\Srp\^ + \Slr\^\Slp\^ + \Slr\^\Srp\'' 



Here, we used the sum rule 



A=L,R,P A=L,R,P 



E 



(9) 



(10) 



which is a direct consequence of the unitarity of the scat- 
tering matrix §. On the other hand, if [S'lpP < |5'pLp, 
we can use (fTO| to write 

^^i\spL\' + \SRLmspL\' + \spRn 

-ISlpI'-ISplI^' + ISlpI^ISplI^ 

= \Spl\^\Spr\^ + \Srl\^\Spl\^ + \SRL\^\SpRf 

+ \Slp\H\Spl\'~\Slp\^)>0. 

(11) 

In conclusion, we have shown that Det V > and there- 
fore the proof that IK' is positive semidefinite is complete. 



